Langevin equation with scale-dependent noise 



M.V.Altaisky 

Joint Institute for Nuclear Research, Dubna, Russia; and 
Space Research Institute RAS, Profsoyuznaya 84/32, Moscow, 117997, Russia 

March 21, 2003 

Abstract 

A new wavelet based technique for the perturbative solution of the Langevin equa- 
tion is proposed. It is shown that for the random force acting in a limited band of scales 
the proposed method directly leads to a finite result with no renormalization required. 
The one-loop contribution to the Kardar-Parisi-Zhang equation Green function for the 
interface growth is calculated as an example. 

The Langevin equation is one of the most general approximations for the evolution of a dy- 
namical system affected by fluctuating environment. It arises in the description of magnetic 
at the presence of magnetic field fluctuations, in the description of hydrodynamic turbu- 
lence, in stochastic quantization problems, in the description of interface growth and in a 
large variety of other problems [TJ 12 El E] ■ 

In the most general form the Langevin equation can be casted as 

?*^ = U[<l>(x,t)]+Ti(x,t) t (r,(x)r](x')) = D(x,x'), (1) 

where U[(f)} is the nonlinear interaction potential, r](t, x) is the Gaussian random noise, 
which standss for the fluctuations of the environment. The Minkovski-like (d+l)-dimensional 
notation x = (x, t),k = (k, u) is used hereafter. 

The standard way to solve the Langevin equation (jj) is to introduce the small parameter 
A in the interaction potential U, and then solve the system iteratively in each order of 
the perturbative expansion. The averaging over the Gaussian random force rj consists in 
evaluation of the pair correlators [rjrj) . The procedure is simplified by the assumption of 
the Gaussian statistics of the random noise which allows one to take into account only even 
order correlators of the random noise: all terms containing the odd number of r\ are equal 
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to zero. The diagram technique for the iterative solution of the Langevin equation is often 
called the Wyld diagram technique |3] . Similarly to that in quantum field theory, it requires 
elimination of loop divergences by renormalization group methods [21 El- 

The structure of divergences arising in the perturbative solution of the Langevin equation 
depends on a particular type of the random force correlation function. Most approximations 
use 5-correlated random force of the form 



The spatial part of the correlation function D(k) is assumed either to be a constant or to 
have a power behavior (D(k) ~ |fc| _/3 ). 

From the physical standpoint, a random force acting in a limited range of scales is often 
desirable. For instance, a limited band forcing is the case for the stirred hydrodynamic 
turbulence and stochastic interface growth. In this paper we propose a novel method to 
describe the limited band stochastic forcing. It is shown, that for a narrow band forcing, an 
appropriate choice of the random force correlation functions yields a theory which is free of 
loop divergences. The proposed method preserves the whole structure of the perturbation 
expansion, with only the space of functions being changed. In the limiting case of the scale 
independent forcing (j2j) all common results are preserved. 

To study the dynamical system separately at each scale, following [6 a , instead of the 
usual space of the random functions f(x,-) G (Q, A, P), where f(x) G L 2 (R d ) for each given 
realization of the random process, we turn to the multi-scale representation provided by the 
continuous wavelet transform 



The wavelet transform is performed here only in the spatial argument, but not in its temporal 
argument of x. This is because the structure of divergences and other important properties 
of the considered processes are determined by their spatial localization. 
The existence an uniqueness of the inverse wavelet transform 



{V(ki)v(h)) = (27i) d+1 5(k 1 + k 2 )D(k 2 ). 



(2) 




(3) 





is provided by the admissibility condition imposed on the basic wavelet 
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Factor Sd, the area of a unit sphere in d dimensions, stands for the rotationally symmetric 
wavelet ip(x) = ip(\x\). 

The use of the wavelet images instead of the original stochastic processes provides an 
extra analytical opportunities. What is important for us, is the possibility of construction 
of more than one set of random functions W(a,b,-), the images of which have coinciding 
correlation functions in the space of f(x, •). 

It is easy to check, that the image of the correlation function of the processes 



Therefore, modeling the random force in the space of wavelet coefficients, we can provide a 
narrow band forcing keeping intact all required properties in the common space. 

As an example, let us consider the well known Kardar-Parisi- Zhang model of interface 
growth 



(W(a u h)W(a 2 , k 2 )) = C42%) d 5 d (k 1 + k 2 )a d+l S( ai - a 2 )D t 



coincides with that of the white noise 



(m)f(k 2 )) 

(Wia^kjWia^h)) 




(27r) d D 5 d (k 1 + fc 2 )(a 1 a 2 ) d / 2 ^(aifc 1 )^(a 2 A; 2 ). 




(6) 



Substituting the wavelet transform in spatial argument 





into the equation (|BjL with the random force of the form 



(77(01, ki)rj(a 2 , k 2 )) 
(fj(a,k)) 



C^nY+^^ih + k 2 )a d+1 5{ ai - a 2 )D(a 2 , k 2 ) 
0. 




after straightforward calculations, we get an integral equation 




ki(k — ki)Z(ai, ki)Z(a 2 , k — ki) 



d d+1 k\ dai da.2 
(2 7 r) d + 1 ai^ 1 a 2 d+1 ' 



from where, the one loop approximation for the Green function is yielded 



G(k) = G (k) - \ 2 G 2 (k) 



d d+1 kx 
(2n) d+1 



A(fci)fci(fc - fc 1 )|G (A; 1 )| 2 fcfc 1 Go(A; - k t ) + 0(A 4 ) (9) 
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where G Q 1 (k) = —iuj + vk 2 is the zeroth order approximation for the Green function, and 

A(k)=C?f ^\ij(ak)\ 2 D(a,k) (10) 

is the scale averaged effective force correlation function. The Green function obtained with 
the random force (JHJ), does not depend on scale Z(a, k) = G(k)fj(a, k). 
Similarly, for the pair correlation function (ZZ), we get 

A 2 — 

C(ai,a f ,k) = —\G {k)\ 2 i)( ai k)i)(-a f k) (11) 

d d+1 h 



/rjd+lh. 
^^ r |Go(fci)| 2 |G (A; - k^Mk - fc 1 )] 2 A(fc 1 )A(A; - h). 



For the random stirring which does not depend on scale, the integration in equation ([lUp. 
after substitution k\ — k' x + | in both equations ©and (fTTj) . leads to the known result [3]. 
Let us consider a single band forcing 

D(a,k) = 5(a-a )D(k) (12) 

and the "Mexican hat" taken as the basic wavelet 

^(jfe) = (27r) d / 2 (-2fc) 2 exp(-fc 2 /2), = (2vr) d . (13) 

Substituting (|T2|) and (fT^j) into the equation 0, after integration over the frequency, in the 
leading order in the small parameter x = |fe|/|fei|, we get the contribution to the Green 
function (d > 2): 

G(k) = G (k) + ^mj^yi^^f I £(<z)e-^Y+Mg + 0(A 4 ). (14) 

For constant D(q) = Dq, the obtained contribution to the Green function is finite and does 
not require any further renormalization. In the limit u, k^O, the one loop contribution to 
the surface tension u, which stems from (114)) . is equal to 



v, 



eff 



X Sd A- d D T l + d -)+0(\* 



v*k 2 16d (27r) d u u V 2 



(15) 



We have to note that the equations ()14I15)) were derived from one loop integral Q after 
substitution ki = q + k/2, followed by perturbation expansion in small parameter x = k/q, 
[7] . For concrete correlation functions the direct numeric evaluation of the integral (jHJ) may 
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be more adequate. Similar calculations can be done for other interaction potentials. For 
instance, for the quadratic interaction \Z 2 , the corresponding equations for (JHJ) and (fTTj) 
differ from that obtained above only by the absence of the scalar products of the wave vectors 
in each vertex. 

Concerning the contribution of the higher orders of the perturbation expansion, we should 
say that for the basic wavelets ip(k) localized in /c-space and limited band noise D(a, k), the 
effective coupling constant A, which is the actual parameter of the perturbation expansion 
jS], can be made small by decreasing the noise amplitude. For instance, for the basic wavelets 
from the family i)(k) = (2 r K) d l 2 (—ik) n exp(— k 2 /2) and the noise correlator -D(a, k) = DoS(a — 
ao), the effective coupling constant is 

\ 2 = ^£± 

z/ 3 a ' 
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